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Abstract 

The Nrf family of transcription factors mediates adaptive responses to stress and longevity, but the identities of the crucial 
Nrf targets, and the tissues in which they function in multicellular organisms to promote survival, are not known. Here, we 
use whole transcriptome RNA sequencing to identify 810 genes whose expression is controlled by the SKN-1/Nrf2 negative 
regulator WDR-23 in the nervous system of Caenorhabditis elegans. Among the genes identified is the synaptic cell adhesion 
molecule n/g-7/neuroligin. We find that the synaptic abundance of NLG-1 protein increases following pharmacological 
treatments that generate oxidative stress or by the genetic activation of skn-1. Increasing nlg-1 dosage correlates with 
increased survival in response to oxidative stress, whereas genetic inactivation of nlg-1 reduces survival and impairs skn-1- 
mediated stress resistance. We identify a canonical SKN-1 binding site in the nlg-1 promoter that binds to SKN-1 in vitro and 
is necessary for SKN-1 and toxin-mediated increases in nlg-1 expression in vivo. Together, our results suggest that SKN-1 
activation in the nervous system can confer protection to organisms in response to stress by directly regulating nlg-1/ 
neuroligin expression. 
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Introduction 

Oxidative stress is generated in cells when an imbalance occurs 
between the production of electrophilic reactive species and the 
endogenous defenses against these harmful molecules [1]. 
Increased reactive oxygen species (ROS) can cause damage to 
cellular components such as lipids, DNA, and proteins, and can 
result in a myriad of detrimental effects, including protein 
aggregation, changes in cell signaling, or altered cell cycle 
progression. Given the nervous system's high metabolic demands, 
high lipid and iron content, and low regenerative ability, it is not 
surprising that oxidative stress is particularly detrimental to this 
tissue; indeed, elevated levels of reactive species have been 
implicated in a range of neurodegenerative diseases including 
Parkinson's disease, Alzheimer's disease and amyotrophic lateral 
sclerosis [2-4] . In multicellular organisms, investigating pathways 
that can mitigate the consequences of elevated oxidative stress on 
a cellular and organismal level is an ongoing area of investigation. 

The Nrf (nuclear factor E2-related factor) family of transcription 
factors controls the primary response to oxidative and xenobiotic 
stress in mammals [5-7] . Studies of Nrf2 have established a critical 
role for the transcription factor in defending tissues from oxidative 
damage [7,8]. During normal conditions, Nrf2 is sequestered in 
the cytoplasm by the kelch-domain containing protein Keapl; in 
response to stress, Keapl releases Nrf2, allowing the transcription 



factor to translocate into the nucleus and initiate transcription of 
downstream targets by binding to the antioxidant response 
element (ARE) in the promoter region of stress response genes 
[6,9]. 

Nrf2 is ubiquitously expressed, and in neurons, increased Nrf2 
activity protects against toxicity by hydrogen peroxide, glutamate, 
and mitochondrial toxins [10,11]. Microarray studies using 
neuronal cultures lacking Nrf2 have identified detoxifying, 
antioxidant and defense genes to be regulated by the transcription 
factor [10-12]. Interestingly, these studies have also found changes 
in the expression of genes involved in a variety of processes, 
including cell signaling and calcium homeostasis, as well as 
neuron-specific genes, but the functional significance of these 
changes has not been determined. 

In C. elegans, the Nrf family homolog SKN-1 [13] confers 
cellular and organismal protection from a variety of environmental 
stressors. Although most studies have examined detoxification and 
stress response in the intestine, a role for SKN-1 in the nervous 
system is emerging. SKN-1 functions in a pair of neurons to 
promote longevity [14], and systemic knockdown of skn-1 increases 
dopaminergic neuron degeneration during methylmercury, alu- 
minum and manganese toxicity [15-17]. Furthermore, loss of skn- 
1 enhances manganese-induced organismal death [18]. In C. 
elegans, the abundance of SKN- 1 is negatively regulated in part by 
the DDB-l/CUL-4 ubiquitin ligase substrate targeting protein 
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Author Summary 

Organisms have evolved mechanisms to protect them- 
selves at the cellular level in response to a variety of 
environmental stresses. Oxidative stress, caused by an 
imbalance in the cellular production of free radicals and 
endogenous antioxidant defenses, can be particularly 
detrimental to the nervous system. Indeed, elevated levels 
of oxidative stress have been linked to nearly all 
neurodegenerative diseases. Therefore, understanding 
how living creatures protect themselves against oxidative 
stress, from a cellular to a systemic level, is vital. Here, we 
have found that increased stress activates a pathway that 
increases the amount of certain proteins found at neuronal 
synapses. This presents an interesting model in which, in 
response to stress, neurons might attempt to enhance the 
strength of the synapse to prevent degeneration. 

WDR-23, which is proposed to function in an analogous manner 
as Keap 1 to promote degradation of SKN- 1 during non-stressed 
conditions. [19,20]. Mutants lacking wdr-23 have increased SKN-1 
protein levels, express high levels of genes involved in antioxidant 
and xenobiotic responses, and are resistant to stress [19-21]. 

In this study, we find that SKN-1 is negatively regulated by 
WDR-23 in cholinergic motor neurons. Using whole transcrip- 
tome RNA sequencing (RNAseq) of wdr-23 mutants expressing 
functional wdr-23 a in the nervous system, we identify 810 genes 
whose expression is likely to be regulated by SKN-1; one of these is 
the cell adhesion molecule nlg-1 /neuroligin. Both nlg-1 expression 
in neurons and NLG- 1 protein abundance at synapses increase in 
mutants with increased SKN-1 activity as well as in animals 
exposed to toxins that generate mitochondrial stress. Furthermore, 
increasing NLG- 1 protein abundance enhances survival following 
toxin treatment, and loss of nlg-1 diminishes SKN- 1 -mediated 
toxin resistance. Together, these results support a role for SKN-1 
in promoting organismal survival by regulating synaptic neuroligin 
abundance. 

Results 

WDR-23 regulates SKN-1 abundance in neurons 

The skn-1 locus encodes three isoforms, skn-la, skn-lb and skn-lc 
(wormbase.org), which differ in their N termini and utilize unique 
transcriptional start sites, skn-lc is primarily detected in the 
intestine; .skn-lb, on the other hand, is expressed principally in a 
pair of sensory neurons and is involved in dietary restriction 
induced longevity [13,14]. The expression of the longest isoform, 
skn-la, however, has not been examined. To examine the 
expression of skn-la, we generated a fluorescent reporter by using 
a 7.3 kb fragment upstream of the skn-la start site to drive gfp fused 
to a nuclear localization signal (nls-gfp). skn-la is a downstream 
gene within the bec-1 /heclin operon [14,22,23], and the 7.3 kb 
fragment includes the bec-1 coding region and 3.3 kb upstream of 
bec-1 (Figure 1 A). In animals expressing this construct, fluorescence 
was observed in many tissues, as previously reported [24,25], as 
well as in cholinergic neurons of the ventral cord (labeled with the 
unc-17/VACKT promoter expressing mCherry) and GABAergic 
neurons (unlabeled ventral cord neurons; Figure IB). 

WDR-23 has two isoforms, WDR-23a and WDR-23b, which 
are expressed in neurons of the ventral cord [19-21]. In 
cholinergic neurons of the ventral cord, a WDR-23a::GFP fusion 
protein was found in cell bodies and in axons where it adopted a 
punctate pattern of fluorescence (Figure 1C), whereas WDR- 
23b::GFP localized exclusively to nuclei (Figure ID). In axons, 



WDR-23a puncta co-localized with a synaptic vesicle marker 
(Figure 1C), and WDR-23a fluorescence became diffuse in animals 
lacking unc-1 04 /kinesin (Figure IE'), which is required for 
trafficking of organelles along microtubules in neurons [26]. This 
suggests that WDR-23a associates with presynaptic organelles. 
WDR-23a also co-localized with a mitochondrial marker at 
synapses (Figure IE), and in unc-104 mutants, WDR-23a puncta 
remained co-localized with displaced mitochondria in axons 
(Figure IE'). In mutants defective in drp-l/T)rp\, a protein 
required for normal mitochondrial fission [27], both WDR-23 and 
mitochondrial markers became displaced in axons (Figure IE"). 
Together, these results suggest that WDR-23a may associate with 
presynaptic organelles, including mitochondria, at synapses. 
Consistent with this, WDR-23a co-localizes with outer mitochon- 
drial membrane markers when expressed in muscle cells [20]. 

To test whether WDR-23 regulates the abundance of SKN-1 in 
the nervous system, as it does in the intestine [19,21], we examined 
the abundance of SKN-la::GFP fusion proteins in motor neurons. 
We drove SKN- la:: GFP expression using the cholinergic unc-171 
VAChT promoter since unc-1 7 is not transcriptionally regulated by 
wdr-23 [20], and therefore changes in SKN-la::GFP should reflect 
changes in protein abundance and not skn-1 expression. In wild 
type animals, SKN-la::GFP fluorescence was detected in 49% 
(n = 272) of ventral cord cholinergic neurons (Figure 2A), consis- 
tent with low levels of SKN-1 reported in the intestine in non- 
stressed animals [13]. In wdr-23 mutants, SKN-la::GFP fluores- 
cence was detected in a larger fraction of neurons compared to 
wild type animals (77%, n = 349). In addition, the average 
fluorescence intensity of SKN- la:: GFP in neuronal cell bodies of 
wdr-23 mutants increased compared to wild type controls 
(Figures 2A and 2B; mean fluorescence ± sem wt: 44.7±2.8; wdr- 
23: 62.6±3.1; /j<0.001, Student's /-test). Together, these results 
suggest that WDR-23 negatively regulates the abundance of SKN- 
la in motor neurons. 

A subset of neuronal genes is regulated by SKN-1 

wdr-23 mutants have a variety of developmental and behavioral 
defects including reduced locomotion, resistance to stress, small 
size, and developmental delay. Each of these defects is suppressed 
by loss of skn-1 [19-21], suggesting that the primary function of 
WDR-23 is to negatively regulate SKN-1. We previously identified 
2,285 transcripts that are significantly up-regulated in wdr-23 
mutants compared to wild type controls using RNAseq [20]. To 
identify genes among these that are regulated by WDR-23 
specifically in the nervous system, we performed RNAseq of wdr- 
23 mutants expressing an integrated array of full-length wdr-23a 
cDNA driven by the pan-neuronal snb-1 promoter (referred to as 
WDR-23 rescue). Of the 2,285 genes whose expression increased 
in wdr-23 mutants, transcripts of 8 1 0 were significantly reduced in 
the WDR-23 rescue animals (Figure 3A and Table SI). The 
average expression levels of the rescued genes in wdr-23 mutants 
was 20.74 fold above wild type which is similar to the average of 
23.99 for all 2,285 genes, indicating that the rescue was not caused 
by a bias generated by rescuing specifically low expressing genes in 
wdr-23 mutants. 

We used the Database for Annotation, Visualization and 
Integrated Discovery (DAVID) [28], to assign GO terms for the 
810 rescued genes; of the 492 genes with predicted functions, 
many genes involved in stress, detoxification and metabolism were 
identified, consistent with the known stress response roles of skn-1 
(Figure 3B). The expression of 177 genes has been experimentally 
examined (wormbase.org), and 93 of these are expressed in the 
nervous system (Figure 3C and Table SI). A subset of these is 
expressed exclusively in neurons (Table 1), including eight genes 
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Figure 1. SKN-1 expression and WDR-23 sub-cellular localization in motor neurons. A) Schematic of bec-1 and skn-1 genetic loci. Locations 
of gfp tags used in this study indicated. B) Co-localization of nls-gfp driven by the 7.3 kb Pbec-1 promoter (vjEx254) with Punc17::rfp (nuls321) in 
cholinergic motor neurons. Arrowheads indicate cholinergic motor neurons; arrow indicates GABAergic neuron. C-D) Left, Localization of WDR- 
23a::GFP {vjEx423; C) and WDR-23b::GFP (vjEx426; D) in ventral cord motor neurons expressed under the snb-1 promoter fragment. Neuronal cell 
bodies are identified by the unc-17 promoter driving rfp {nuls321) for reference. Right, WDR-23a::GFP [vjls26; C) and WDR-23b::GFP (vjEx7; D) co- 
localization with the synaptic vesicle associated protein mCherry::SNB-1 [vjEx339] in motor neuron axons using the unc-129 promoter fragment. E) 
WDR-23a::GFP {vjls26) co-localization with the mitochondrial marker (INVOIvkRFP; vjEx663) in motor neuron axons in wild type, unc-704Ainesin 
mutants (E') and drp- 7/Drp1 mutants (E"). Scale bars represent 10 |im. 
doi:1 0.1 371 /journal.pgen.1 0041 00.g001 



encoding insulins or FMRFamide related peptides, as well as three 
genes encoding proteins involved in neuropeptide processing — egl- 
5/proprotein convertase, egl-21 /carboxypeptidase, and sbt-2/7B2. 
We also identified five genes encoding synaptic proteins, including 
unc-13 /Muncl3, nig- 1 /neuroligin, pde-4/cyclic nucleotide phos- 
phodiesterase, dlk-1 /MAPKKK, and raZ>-7/NPDC-l [29-32]. 
Together, these results suggest that neuronal WDR-23, possibly 
through SKN-1, regulates the expression of genes involved in 
neuroendocrine signaling and synaptic function. 

To identify genes that may be direct SKN- 1 targets, we examined 
the promoters of the neuronal genes for the presence of consensus 
SKN-1 binding sites. SKN-1 is predicted to bind the consensus 
sequence WWT DTCAT on either strand in the promoters of target 
genes. Using the web based program Regulatory Sequence Analysis 
Tools (RSAT), we scanned 1000 bp promoter fragments for each 
gene for the SKN-1 consensus [33]; the promoter fragments for 70 
of the 93 neuronal genes contained at least one potential binding 
site (Table SI). We then cross-referenced the list of rescued 



neuronal genes with SKN-1 ChlP-seq datasets taken from LI, L3 
and L4 stage animals (ModEncode and [34]) and found that 31 
genes contained significant SKN-1 peaks within 2000 bp upstream 
or 500 bp downstream of the transcription start site (Table SI). 
These results indicate that a subset of the neuronal genes identified 
by RNAseq may be direct binding targets of SKN- 1 . 

nlg-1 expression is positively regulated by SKN-1 

Among the neuronal genes containing SKN- 1 consensus sites in 
their promoter, the cell adhesion molecule nlg-1 /neuroligin 
emerged as an interesting candidate since it has been implicated 
in C. elegans stress response and in synaptic function [30,35,36]. nlg- 
1 is the sole neuroligin family ortholog in C. elegans, and NLG- 1 is 
expressed in cholinergic motor neurons where it localizes to 
presynaptic terminals [35,37]. 

To examine the effects of skn-1 on nlg-1 expression, we 
constructed reporter strains consisting of a 3.6 kb nlg-1 rescuing 
promoter fragment [30] driving soluble gfp (Pnlg-lr.gfp, vjls47, and 
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Figure 2. WDR-23 regulates SKN-1 abundance in motor neurons. A) Representative images of SKN-1a::GFP driven by the unc-17 promoter 
(vjls45) in cholinergic motor neurons of wild type, left, and wdr-23 mutants, right. Cell bodies identified by presence of RFP expressed under the unc- 
17 promoter (nuls321). B) Quantification of SKN-1 a::GFP (vjls45) in wild type (N = 28) and wdr-23 (N = 33) mutants. Cells were classified as low, medium 
or high expressing neurons based on average SKN-1 a::GFP fluorescence per cell (see methods). Scale bar represents 10 urn. 
doi:1 0.1 371 /journal.pgen.1 0041 00.g002 



vjls48; Figure 4A). In wild type animals, we detected Pnlg-lr.gfp 
fluorescence in a few head neurons and in ventral nerve cord 
neurons, previously reported to be DA and VA class cholinergic 
motor neurons [30,35]. nlg-1 reporter fluorescence has been 



reported at low levels in muscles [30], but we did not detect 
Pnlg-l::gfp in muscle cells in either transgenic line. In mutants 
lacking wdr-23, Pnlg-lr.gfp fluorescence in motor neurons increased 
approximately 5 fold (Figure 4B and 4C), in agreement with the 
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Figure 3. Analysis of whole transcriptome RNA sequencing of wdr-23 mutants expressing neuronal WDR-23. A) Correlation plot of 
RNAseq reads per kilobase per million reads (RPKM) for WDR-23 rescue compared to wdr-23 mutants, for 2,285 genes previously identified to be 
significantly different between wild type and wdr-23 mutants [20]. Correlation plot includes 2,1 18 genes, as successful reads were not obtained for all 
genes in the rescuing line. WDR-23 rescue denotes wdr-23 mutants with expression of WDR-23a (using snb- //synaptobrevin promoter, nuls225). 
Genes that do not rescue (indicated in blue) have similar RPKM for wdr-23 mutants and the WDR-23 rescue, whereas genes that rescue (indicated in 
red) have statistically higher RPKM for wdr-23 mutants. B-C) Characterization of 810 genes significantly different in wdr-23 versus WDR-23 rescue 
(nuls225). B) Right, Breakdown by GO term of genes with predicted functions identified by DAVID (492 genes). Left, Breakdown by GO term of rescued 
genes expressed specifically in the nervous system (43 genes). C) Distribution by tissue expression of genes with WormMart reported expression 
patterns (1 73 genes). 
doi:1 0.1 371 /journal.pgen.1 0041 00.g003 
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Table 1. Neuronal targets of SKN-1 identified by RNA sequencing. 
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RNAseq list of genes with neuronal expression patterns according to WormMart WS220 that were rescued by Psnb-1::wdr-23a in wdr-23 mutants. Functions identified 
by DAVID. Consensus indicates the number of SKN-1 binding sites (WWTDTCAT) in promoter fragments within 1 kb of transcriptional start site as evaluated by RSAT. 
doi:1 0.1 371 /journal.pgen.1 0041 OO.tOOl 



5.6 fold increase in nlg-1 transcript levels detected by RNAseq 
[20]. Cell-specific expression of either wdr-23a or wdr-23b cDNA in 
nlg-1 expressing cells (using the nlg-1 promoter) fully rescued the 
increased Pnlg-l::gfp fluorescence to wild type levels (Figures 4B 



and 4C), indicating that the regulation of nlg-1 expression by wdr- 
23 is cell autonomous. 

In skn-1 mutants lacking the skn-la/c isoforms (zu67 mutants), 
Pnlg-l::gfp expression was similar to wild type. However, skn-1 is 
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Figure 4. nlg-1 transcription is regulated by neuronal SKN-1. A) Schematic of Pnlg-1::gfp reporters; vjls47 and vjls48 contain a 3.6 kb nlg-1 
promoter driving gfp. Location of SKN-1 binding site and mutated site indicated. B) Top, Representative images of Pnlg-1::gfp {vjls47) expression in the 
ventral cord of indicated strains. WDR-23a and WDR-23b rescues indicate isoform specific cDNA driven by the nlg-1 promoter {vjEx447 and vjEx436, 
respectively) and expressed in wdr-23 mutants. Bottom, Representative images of the extra-chromosomal transgene Pnlg-1 (Abs)::gfp {vjEx391) in the 
indicated strains. C) Quantification of Pnlg-1 ::gfp (vjls47) and Pnlg-1(Abs)::gfp (vjEx391) in indicated strains. Sample sizes indicated. Values normalized 
to wild type (vjls47) or deleted binding site (vjEx391), respectively. D) Representative images of Pnlg-1 ::gfp (vjls48) heterozygotes in indicated strains. 
skn-1 indicates zu67 loss of function allele, unless otherwise indicated. E-F) Normalized quantification of Pnlg-1::gfp (vjls48) heterozygotes (E), Pnlg- 
1::gfp (vjls48) homozygotes and Pnlg-1 (Abs)::gfp (vjEx756, F) in indicated strains. All nlg-1 reporters were imaged using the same microscope settings, 
and expression of reporters was similar between strains (total fluorescence vjls48 in wild type animals: 3280.6± 155.6; total fluorescence of vjEx756 in 
wild type animals: 2919.7±241.7). G) Electrophoretic mobility shift assay for binding of full-length SKN-1 a to the consensus sequence at -396 bp. 
Lysate was either unprogrammed control or programmed to express SKN-1 a. Competitor is 20, 100 and 200 fold molar excess of unlabeled probe. 
Scale bar represents 10 |im; error bars represent ±sem; ***p<0.001. 
doi:1 0.1 371 /journal.pgen.1 0041 00.g004 



required for the increased Pnlg-1 r.gfp fluorescence caused by loss of 
wdr-23, since skn-1 mutations reduced the Pnlg-l::gfp reporter 
fluorescence of wdr-23 mutants to wild type levels (Figures 4D 
and 4E). Conversely, Pnlg-l::gfp fluorescence increased by ~35% 
in mutants in which skn-1 is hyperactive (laxl20gf or laxlSSgf, 
Figure 4F). Iaxl20gfnnd laxlSSgf 'are thought to prevent SKN-la/c 
interaction with mitochondrial docking proteins, resulting in 
an activated pool of SKN-la/c [38]. Together, these results 
show that skn-1 is not necessary for baseline nlg-1 expression in 
motor neurons, but skn-1 activation positively regulates nlg-1 
expression. 

The nlg-1 promoter has four SKN-1 binding consensus sites within 
500 bp upstream of the transcriptional start site (Figure SI). Promoter 
alignments between the nematode Caenorhabditis species elegans, 
briggsae, japonica and remand revealed that one of these sites located 
396 bp upstream of the C. elegans start is completely conserved in all 
species (Figure SI). This site has the sequence AAT GTCAT , which 
matches the consensus perfectly. The underlined region is predicted 
to be a largely invariant sequence that directly interacts with SKN- 1 
[13,39]. We mutated AAT GTCAT at -396 to AACTGCAG in the 
Pnlg-1 r.gfp reporter to create a reporter with a deleted binding site (the 
Pnlg-1 (Absjrgfp reporter; Figure 4A). Basal motor neuron fluorescence 
of transgenic animals expressing Pnlg-1 (Absjrgfp was similar to 
transgenic animals expressing the Pnlg-l::gfp reporter (Figure 4B). 
However, Pnlg-1 (Absjrgfp reporter fluorescence did not increase in 
either wdr-23 mutants or in skn-1 (gf) mutants compared to wild type 
controls (Figure 4B, 4C and 4F), suggesting that this site is critical for 
-mediated increases in nlg-1 expression. To test whether SKN-1 
binds to this site, we performed electrophoretic mobility shift assays. 
In vitro translated SKN- la bound to labeled probes containing the 
putative SKN-1 binding site in the nlg-1 promoter, and binding was 
disrupted by the addition of excess unlabeled probe (Figure 4G). 
These results indicate that the SKN-1 binding site at —396 in the nlg- 
1 promoter can be bound by SKN-1 in vitro and is critical for skn-1- 
induced expression of nlg-1 in vivo. 

Synaptic abundance of NLG-1/neuroligin levels is 
regulated by SKN-1 signaling 

To determine whether the transcriptional regulation of nlg-1 by 
SKN-1 impacts NLG-1 protein levels at synapses, we examined 
synaptic levels of NLG- 1 in animals expressing a fusion protein in 
which GFP was inserted near the C-terminus of NLG-1 (NLG-1 - 
GFP, vjEx561 and vjlslOS). This fusion protein is functional [30] 
and localizes to presynaptic terminals in motor neurons [35,37]. 
NLG-1 -GFP driven by the nlg-1 promoter adopted a punctate 
pattern of fluorescence in the dorsal and ventral cords, where 
presynaptic terminals of DA and VA class motor neurons are 
located, respectively (Figure 5A). We examined changes in the 
synaptic abundance of NLG-1 -GFP by measuring the average 
punctal fluorescence intensity (peak fluorescence) and synapse 
number (interpunctal interval) [40,41]. In skn-1 (gf) mutants, the 



punctal fluorescence of NLG-1 -GFP significantly increased in both 
the dorsal and ventral cords, while the interpunctal interval did not 
change (Figure 5A, 5B and Table S2). These results indicate that 
SKN-1 positively regulates synaptic NLG-1 protein abundance 
but does not affect synapse number. 

Environmental stress induces NLG-1 expression 

We next tested whether toxins that induce oxidative stress could 
increase neuroligin expression and abundance in motor neurons. 
We found that exposure to the mitochondrial stressors juglone or 
sodium arsenite, both of which have been shown to activate skn-1 
[42,43], robustly induced Pnlg-lr.gfp fluorescence in motor neurons 
compared to untreated animals (Figure 6 A and 6B). Treatment 
with an organic mercury (thimerosal) increased nuclear SKN- 
1::GFP in the intestine (Figure S2) and also increased Pnlg-1 r.gfp 
fluorescence (Figure 6B). However, these toxins had no effect on 
Pnlg-1 (Absjrgfp fluorescence (Figure 6A and 6B). We found that 
juglone treatment significantly increased punctal fluorescence of 
NLG-1 -GFP, without changing the interpunctal interval 
(Figure 6C, 6D and Table S2). These results suggest that 
activation of skn-1 by oxidative stress increases NLG-1 synaptic 
abundance by increasing nlg-1 expression in neurons. 

Genetic regulation of nlg-1 expression 

To further explore how SKN-1 activity is regulated in the 
nervous system, we tested the impact of altering insulin signaling, 
mitochondrial respiration or synaptic activity on Pnlg-lr.gfp 
expression. The insulin-like signaling (IIS) pathway regulates 
SKN-1 in the intestine [44]; activation of DAF-2 /insulin-like 
receptor leads to SKN-1 phosphorylation by SGK-1/SGK, 
resulting in decreased SKN- 1 activity. We examined putative null 
sgk-1 mutants and found that Pnlg-lr.gfp reporter fluorescence 
increased by ~25% (Figure S3). Conversely, during conditions of 
stress, PMK-l/p38 MAPK phosphorylates SKN-1, resulting in 
nuclear SKN-1 translocation [45]. We found no change of 
baseline Pnlg-lr.gfp fluorescence in mutants lacking either sek-1 1 
MAPKK or pmk-1 (Figure 6E, 6F and S3). However, sek-1 and 
pmk-1 mutations suppressed the juglone-induced increase in Pnlg- 
lr.gfp reporter fluorescence (Figure 6E and 6F). 

C. elegans mutants with impaired mitochondrial respiration, for 
example the conserved clk-1 /COQ7 , necessary for the biosynthesis 
of coenzyme Q, and isp-l/l&Y, a subunit of mitochorial complex 
III, are resistant to toxins that increase oxidative damage [46] . pink- 
1 /PINK 1 , on the other hand, is predicted to act in conjunction with 
the E3 ligase Parkin to initiate mitophagy of damaged mitochondria 
[47,48]. Pnlg-lr.gfp expression increased mildly in mutants lacking 
clk-1, but not isp-1. Conversely, loss of pink-1 resulted in a significant 
decrease of Pnlg-lr.gfp fluorescence (Figure S3). 

Finally, in order to test whether neuronal activity regulates SKN- 1 
in the nervous system, we examined mutants with increased synaptic 
activity (dgk-1 /diacylglycerol kinase or goa-1 /Goto) and mutants with 



PLOS Genetics | www.plosgenetics.org 



7 



January 2014 | Volume 10 | Issue 1 | e1004100 



Stress Regulation of Neuroligin 




Figure 5. Synaptic NLG-1 levels are regulated by SKN-1. A) Representative images of Pnlg-1::nlg-1-gfp (vjlslOS) in the dorsal (fop) and ventral 
(bottom) cords of the indicated strains. skn-1(gf) indicates Iax188 gain-of-function allele. B) Quantification of punctal fluorescence (peak), cord 
fluorescence (cord) and interpunctal interval (IPI) of Pnlg-1::nlg-1-gfp (vjls105) in dorsal (wild type n = 31, skn-1(gf) n = 29) and ventral cords (wild type 
n = 32, skn-1(gf) n = 26) of the indicated strains. Cord fluorescence could not be measured for vjls105 in the dorsal cord because it was not above 
background levels. Scale bar represents 10 urn; error bars represent ±sem; ***p<0.001. 
doi:1 0.1 371 /journal.pgen.1 0041 00.g005 



decreased activity (unc-2 /VGCC or unc- 1 8 /nSecl), as well as mutants 
lacking mef-2/MEF and mtV-7/microRNA, which are involved in a 
retrograde synaptic signaling pathways that is dependent on nlg-1 
[35]. We detected no change in reporter fluorescence in these 
mutants (Figure S3). Together, these results indicate that insulin 
signaling and mitochondrial metabolism contribute to the activation 
of SKN- 1 in neurons, while changes in synaptic transmission do not 
seem to impact neuronal SKN-1 activity. 

NLG-1 promotes organismal survival during SKN-1- 
dependent response to mitochondrial stress 

We next sought to determine whether nlg-1 mediates the protective 
effects of skn-1 activation in response to environmental toxins, skn-1 
mutants are hypersensitive to toxicity induced by arsenite and juglone 
treatment (Figure 7A and S4; [42,43]). In addition, we found that skn- 
1 mutants were sensitive to thimerosal-induced toxicity (Figure S5), in 
agreement with studies showing SKN- 1 protects from metal toxicity 
[15,18]. In contrast, animals lacking wdr-23 were resistant to toxicity 
of both thimerosal and juglone, and resistance was completely 
blocked by loss of skn-1 (Figures 7A and S5; [43]). Similarly, 
hyperactive skn-1 (gf) mutants were resistant to juglone toxicity 
(Figure 7B and S6). Interestingly, loss of wdr-23 did not confer 
protection against sodium arsenite (Figure S4), suggesting specificity 
in drug resistance in mutants lacking wdr-23. 

nlg-1 mutants are hypersensitive to heavy metal toxicity by 
thimerosal and oxidative stress induced by paraquat [30]. We 
found that mutants lacking nlg-1 were also more sensitive to 
juglone toxicity (Figure 7A and S5). In contrast, transgenic animals 
over-expressing NLG-1 -GFP were significantly more resistant to 
juglone-induced toxicity compared to non-transgenic controls 
(Figure 7C). To confirm that the juglone resistance caused by 
NLG-1 -GFP transgenes was due to nlg-1 expression, we examined 
juglone responses of animals expressing Pnlg-l::gfp and found that 
they were not as resistant to juglone as NLG-1 -GFP expressing 
animals (Figure S7). Finally, nlg-1 mutations dramatically reduced 
the ability of activated skn-1 to protect animals from the toxic 
effects of juglone: skn-1 (gf);nlg-l double mutants were significantly 
less resistant to juglone than skn-1 (gf) mutants alone (Figure 7B). 
Together, these data suggest that the dosage of nlg-1 is a critical 



determinant of survival in response to stress, and that nlg-1 
contributes to .^-/-dependent survival. 

Discussion 

In this study, we find that WDR-23 in neurons, possibly acting 
through SKN-1, regulates the expression of approximately 800 
genes in multiple tissues, including the nervous system. One of 
these genes, the synaptic cell adhesion molecule nlg-1 /neuroligin, 
contains a consensus SKN- 1 binding site in its promoter necessary 
for SKN-1 binding in vitro and for nlg-1 expression following SKN- 
1 activation. Increasing SKN-1 activity elevates the abundance of 
NLG-1 protein at synaptic terminals, and changes in synaptic 
NLG-1 levels correlate with altered survival of animals in the 
presence of mitochondrial toxins. We propose a model in which, in 
response to mitochondrial stress, SKN-1 induces transcription of 
nlg-1, which in turn enhances organismal survival (Figure 7D). 

SKN-1 transcriptional targets in neurons 

Previous studies have demonstrated SKN-1/Nrf2 dependent 
transcriptional programs are initiated in response to oxidative and 
xenobiotic stress, and these programs are critical for organismal 
survival and longevity [42,43,49]. Among the genes regulated by 
SKN-1, few have been shown to mediate the protective effects of 
SKN-1, and fewer still have been shown to be direct binding 
targets. Furthermore, attempts to identify a comprehensive set of 
genes regulated by SKN- 1 /Nrf2 required to protect organisms from 
stress have been difficult due to the lack of tissue-level resolution. 

Here, we have used comparative whole transcriptome RNA 
sequencing to identify an inclusive set of genes that are likely to be 
regulated by SKN-1 in the nervous system. For this analysis, we 
examined wdr-23 mutants expressing a rescuing wdr-23a transgene 
driven by the pan-neuronal snb-1 promoter. Genetic studies 
indicate that several distinct phenotypes displayed by wdr-23 
mutants are completely suppressed by skn-1, indicating that SKN- 
1 is selectively activated in wdr-23 mutants, and a yeast-two hybrid 
screen identifies the only binding target of WDR-23 to be SKN-1 
[19], indicating that SKN-1 is selectively activated in wdr-23 
mutants. Interestingly, wdr-23 mutants appear to activate SKN-1 
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Figure 6. nig- 1 transcription is regulated by the SKN-1 pathway. A-B) Representative images (A) and quantification (B) of animals expressing 
the Pnlgl::gfp {vjls47) or Pnlg-1 (Abs)::gfp (v]Ex756) reporters exposed to control, 120 |iM juglone, 5.0 mM sodium arsenite or 35 uM thimerosal. L4 
stage animals were exposed to drug overnight for 14 hours, followed by 2-4 hour recovery prior to imaging. Values of stressed animals were 
compared to non-stressed animals of the same genotype to determine significant changes in response to stress. Sample sizes indicated. C-D) 
Representative images (C) and quantification (D) of Pnlg-1::nlg-1-gfp (vjls105) in the dorsal cord of animals exposed to control (n = 30) or 120 |j.M 
juglone (n = 29) treatments. E-F) Representative images (E) and quantification (F) of Pnlg-1 ::gfp (vjls47) in indicated strains exposed to control or 
120 uM juglone. Sample sizes indicated. Scale bar represents 10 |im; error bars represent ±sem; **p<0.01, ***p<0.001. 
doi:1 0.1 371 /journal.pgen.1 0041 00.g006 



to a greater extent than skn-l(gf) mutations or toxin treatment. For 
example, loss of wdr-23 results in greater survival in response to 
juglone treatment than skn-l(gf) mutations. In addition, wdr-23 



mutations increase nlg-1 expression in the ventral cord neurons 
fivefold, whereas skn-l(gf) or toxin treatment increased it by 
approximately 30%. Thus, wdr-23 mutants may provide increased 
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Figure 7. NLG-1 is protective to juglone. A-B) Survival curves of indicated strains on 200 uM juglone. C) Survival curve of indicated strains on 
400 uM juglone. D) Model for SKN-1 -dependent transcription of nlg-1. Stress activated the SKN-1 pathway, which increases nlg-1 expression; 
increased NLG-1 activity promotes organism survival. Error bars represent ±sem. Log rank tests were used to identify significant differences between 
genotypes and reported in Table S3. Curves depict average values for four replicates of n = 40 per genotype. 
doi:1 0.1 371 /journal.pgen.1 0041 00.g007 



sensitivity when used for transcriptional profiling, maximizing our 
ability to identify SKN- 1 targets expressed in low abundance or in 
a small subset of cells. It is, however, possible that WDR-23 has 
functions beyond SKN-1 regulation, in which case some of the 
genes identified by this approach may not be targets of SKN- 1 . 
While the snb-1 promoter fragment drives expression of GFP 
strongly in the nervous system, the possibility that this promoter 
may be leaky raises the prospect that some of the genes identified 
here may be regulated by WDR-23 in other tissues in addition to 
the nervous system. 

The 810 genes identified here most likely represent either direct 
SKN-1 targets or indirect targets of SKN-1 that are secondarily 
activated in neurons or other tissues. SKN-1 ChlP-seq of larval 
stage animals identified a list of approximately 3000 genomic 
peaks bound by SKN-1 in vivo [34]; we cross-referenced our gene 
list of rescued neuronal genes with SKN-1 ChlP-seq datasets taken 



from LI, L3 and L4 stage animals [34] and found that a subset of 
the genes we identified contain a significant SKN-1 peak near to 
their transcriptional start sites (Table SI), but many genes do not. 
Differences between these datasets may be a byproduct of using 
non-stressed animals for the ChlP-seq experiments, reflecting 
basal, but not stress-induced, promoter occupancy by SKN- 1 . 

Among the genes we identified were several neuropeptides and 
insulins, as well as peptide processing enzymes reported to be 
expressed in neurons. We speculate that stress-induced peptide 
processing and release may be part of a humoral response to 
promote organismal survival. Identification of the precise peptid- 
ergic signaling pathways will help to elucidate the mechanisms by 
which SKN-1 confers survival. We also identified a handful of 
known synaptic genes, including two cell adhesion molecules, nlg-1 
and ncam-1, indicating SKN-1 might play a role in maintaining the 
stability of neural networks. A theory has emerged suggesting 
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axonal retraction might precede neuronal apoptosis in neurode- 
generation, called dying back degeneration [50], and synaptic 
breakdown may occur prior to axonal retraction. The identifica- 
tion of cell adhesion molecules in this study suggests that prior to 
synaptic breakdown, neurons might initiate transcriptional pro- 
grams to protect the integrity of the synapse. In support of this, 
recent evidence suggests that increased expression of NCAM in 
SH-SY5Y cultures prevents oxidative stress-induced apoptosis, 
and over-expression of a truncated NCAM molecule protects 
neuronal tissue in lesioned rats [51,52]. 

SKN-1 as a mitochondrial stress sensor 

Activation of SKN-1 /Nrf2 can protect against cell toxicity 
induced by oxidative stress [53-55]. In primary neuronal cultures, 
for example, activation of Nrf2 prevents cell death in response to 
rotenone and MPP, potent inhibitors of mitochondrial respiration 
[11]. Loss of OPA1, a key regulator of the morphology of 
mitochondria, results in increased Nrf2 activation [56]. These 
studies and others collectively suggest a potential role for SKN-1/ 
Nrf2 as a sensor for increased mitochondrial dysfunction in the 
nervous system. We found that genes involved in drug detoxifi- 
cation, including the glutathione precursors gst-4, gst-10, gcs-1, and 
gst-1 were reduced by neuronal expression of wdr-23a. Interest- 
ingly, gst-1 has also been shown to contribute to dopaminergic 
neuron survival after manganese treatment [16]. Given their 
confirmed expression in neurons, is possible that these genes have 
a role in direcdy enhancing neuronal protection. 

SKN- 1 associates with purified mitochondrial fractions, and the 
skn-l(gf) alleles are proposed to reduce mitochondrial association, 
suggesting that mitochondria may act as a sink for SKN-1 [38]. 
Our data is consistent with the idea that WDR-23a associates with 
presynaptic organelles, including mitochondria. First, WDR-23a 
localizes to presynaptic terminals, where mitochondria are 
abundant. Second, WDR-23a remains associated with mitochon- 
dria in mutants in which mitochondria are displaced. Third, 
WDR-23a localizes to the outer membrane of mitochondria in 
muscle cells [20]. We speculate that WDR-23 may be a 
mitochondrial stress sensor that regulates SKN-1 abundance. 

Neuroligin as a neuronal target of SKN-1 

Our data supports the idea that SKN-1 activation leads to 
increased NLG-1 abundance at synapses. We found that nlg-1 
confers some, but not all, of the protective effects of activation of 
SKN-1, as nlg-1 mutations reduced, but did not eliminate, the 
resistance of skn-l(gf) mutants to juglone. Additional SKN-1 
targets either in neurons or in other tissues are likely to contribute 
to organismal survival in response to juglone treatment. Interest- 
ingly, loss of nlg-1 did not suppress wdr-23 mutants resistance to 
thimerosal or juglone (Figure S5 and data not shown). This may be 
due to higher SKN- 1 activity in wdr-23 mutants compared to skn- 
l(gf) mutants, which may compensate for the lack oinlg-1. Because 
nlg-1 mutants themselves are hypersensitive to stress, it is possible 
that nlg-1 and skn-1 function in parallel pathways to promote 
resistance; however, our data do not support this idea, but rather 
support the notion that nlg-1 is a direct transcriptional target of 
SKN-1. Basal nlg-1 transcription is not likely to be under skn-1 
regulation, since nlg-1 reporter expression remains unchanged in 
skn-1 mutants. Consistent with this, SKN-1 does not occupy the 
SKN-1 binding site we identified at position —396 in the nlg-1 
promoter in unstressed larval animals by ChlP-seq analysis. 

How might increased expression of a synaptic cell adhesion 
molecule promote organismal survival in response to stress? 
Recent work has established the presence of "mitokines" in 
neurons — a signal produced in the neurons in response to 



dysfunctional mitochondrial electron transport [57]; release of 
these mitokines results in increased organismal survival. It is 
possible that nlg-1 may be required in neurons for proper release of 
mitokines. Furthermore, the neuropeptides identified in this study 
are candidates for being signals released in response to stress to 
promote resistance, nlg-1 expression has been detected in head 
neurons, motor neurons, and muscle cells [30]; thus, it is possible 
that nlg-1 functions in any of these tissues to convey protection. 

In mammals, neuroligin is a post-synaptic cell adhesion 
molecule that binds to the presynaptic protein neurexin; this 
junction is necessary for maintaining mature synaptic connections 
and normal synaptic transmission [58,59]. In humans, rare 
mutations in neuroligin are associated with autism and other 
cognitive disorders. Some of these mutations reduce neuroligin 
delivery to the cell surface, interfering with synapse development 
and synaptic transmission [60-62]. Interestingly, it has been 
suggested that oxidative stress and mitochondrial dysfunction may 
play a role in the pathogenesis of autism, as certain biomarkers for 
oxidative stress are elevated in autistic patients [63-65]. It is 
possible that individuals with these mutations are unable to 
increase synaptic neuroligin levels in response to stress, and it will 
be interesting to identify the cellular and molecular mechanisms 
underlying neuroligin-dependent survival in response to stress. 

Materials and Methods 

C. elegans strains 

Strains were cultured at 20° using standard methods. All 
experiments were performed on young adult hermaphrodites 
unless otherwise indicated. The following strains were provided by 
the Caenorhabditis Genetics Center, which is funded by the NIH 
National Center for Research Resources (NCRR): sek-l(km4), pmk- 
l(km25), sgk-l(ok538), skn-l(zu67), nlg-1 (ok259), clk-l(e2519), mef- 
2(gvl), mir-l(gk276), isp-1 (gml50), and pink-l(ok3538). Strain wdr- 
23 (tml 817) was provided by the National BioResource Project 
(Japan). The wild type reference strain was N2 Bristol. The 
following strains were also used: drp-1 (tml 108), dgk-l(nu62), goa- 
l(sa734), mc-2(ljl), unc-18(md299), glo-1 (zu391), unc-104(el265), 
skn-1 (laxl20gf), skn-1 (laxl88gf), nuhl52[Pttx-3::RFP, Punc- 
129::GFP::SM1] II, nuls321 [Pmyo-2::GFP, Punc-1 7::mCherry\ , 
nuIs225\Pmyo-2::GFP, Psnb-l::WDR-23a], idls7\rol-6(sul006), Pskn- 
l::skn-lb/c::GFP\, yuIs25[Pmyo-2::GFP, Punc-129::mito-GFP\V, and 
vjIs26[Pmyo-2::GFP, Punc-1 29::WDR-23a::GFP\lll. Mutant strains 
were outcrossed a minimum of 4 times; all integrants were 
outcrossed at least 8 times. 

Molecular biology 

C. elegans cDNA was used to clone all genes into pPD49.26 using 
standard molecular biology techniques, unless otherwise noted. 
Promoter elements were amplified from mixed stage genomic 
DNA. The following plasmids were generated: pTS147[/ J «nc- 
129::invom::rfp], pTS3\[Pbec-l:mls-gfp], pDS284[Pnlg-l::gfp], pDS286 
[Pnlg-l(ABS)::gfp], pKG8[Pnlg-l:mlg-l-gfp], pTS79[Punc-17::skn- 
la::gfp], pDS 1 39 [Punc-i 29::snb- 1 ::mCheny\, pT)S237 [Psnb-l::wdr- 
23a::gfp] , pTS85 [Psnb-l:wodr-23b::gfp\ , pDS334 [Pnlg-l::wdr-23a] , 
pGS335[Pnlg-l::wdr-23b], and pTS199[T7::skn-la]. 

Oligo sequences: 

Pbec-1 

oTS44: ccccccGGATCCcgacaattatacatgttcccc 
oTS45: ccccccGCTAGCcgactgactggattatgatagatcc 
Pnlg-1 

oDS694: ccccccGCATGCtaagcccccgtacgctaacacc 
oXL13: ccccccGGATCCgcctgttcacttccaaattcgc 
Pnlg-1 (Abs) 
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oDS692: cctgttgccccccaaatgCTGCAGttacctcttttcctcccttctacc 
oDS693 ggtagaagggaggaaaagaggtaaCTGCAGcatttggggggcaa- 
cagg 

Transgenic lines 

Transgenic strains were generated by injecting N2 with 
expression constructs (2.5-90 ng/uL) and the co-injection marker 
KP#708 {Pttx-3::rfp, 40 ng/ul) or KP#1106 (Pmyo-2::gfp, 10 ng/ 
|XL). Microinjection was performed using standard techniques as 
previously described [66]. At least three lines for each transgene 
were examined for expression, and representative lines were 
quantified. The following strains were made: vjEx7[Punc-l 29::wdr- 
23b::gfp] vjEx663\Punc-129::invom::rfp\, vjEx254[Pbec-l ::nls-gfp], 
vjEx391 [Pnlg-l(Abs)::gfp] , vjEx756[Pnlg-l (Abs)::gfp] , vjEx561 [Pnlg- 
1 '-'-nlg-1 -gfp\ , vjEx339\Punc-129::snb-l:ynCherry\ , vjEx423\Psnb- 1 ::wdr- 
23a::gfp] , vjEx426[Psnb-l::wdr-23b::gfp] , vjEx447[Pnlg-l::wdr-23a\ , 
vjEx436[Pnlg-l::wdr-23b], vjI.s45[Punc-17::skn-la::gfp\ll, vjIs47[Pnlg- 
l::gfpW, vjE48[Pnlg-l::gfp]l, vjM05\Pnlg-l::nlg-l-gp\lll- 

Microscopy and analysis 

To image animals, adult worms were paralyzed using 2,3- 
butanedione monoxime (BDM, 30 u.g/uL; Sigma) and mounted 
on 2% agarose pads for imaging. Images were captured with a 
Nikon eclipse 90i microscope equipped with a Nikon PlanApo 
60 x or lOOx objective (NA= 1.4) and a PhotometricsGoolsnap 
ES2 camera. For fluorescence imaging of synapses, images were 
captured either from the ventral or dorsal cord near the posterior 
gonadal bend of the worm, as indicated. Images of animals 
expressing nlg-1 reporters were captured at the ventral cord near 
the posterior gonadal bend of the worm. We found that skn-1 
mutants expressing fluorescent integrants arrest at larval stages 
when the integrant is homozygous, making quantification of 
fluorescence markers challenging; however, skn-1 mutants express- 
ing heterozygous integrants develop fully into adults. 

Metamorph 7.0 software (Universal Imaging/Molecular Devices) 
was used to capture serial image stacks, and the maximum intensity 
projection was used for analysis of the dorsal and ventral cords. Line 
scans of the maximum intensity projection image were also 
recorded using Metamorph. The fluorescence intensity values were 
then quantified using Puncta 6.0 software written with Igor Pro 
(Wavemetrics), as previously described [67]. For all experiments, 
fluorescence values were normalized to the values of 0.5 urn 
FluoSphere beads (Invitrogen) captured during each imaging 
session. This was performed to provide a standard for comparing 
absolute fluorescence levels between animals from different sessions. 

To quantify changes in neuronal SKN-la::GFP, an anterior and 
posterior image was taken for each animal in both RFP and GFP 
channels. Cell bodies were identified by the presence of soluble 
mCherry. Average cell body fluorescence was calculated by 
outlining the entire cell body and taking the average intensity. 
Background values were determined by finding the average 
fluorescence of the area immediately adjacent to the cell body and 
were subtracted from each cells' average fluorescence. Cells were 
categorized as low, medium, or high expressing cells using 
arbitrary cut off levels after background subtraction — low 
expressing cells were those cells below the level of detection (less 
than 30 units different than background fluorescence), medium 
expressing cells were between 30-70 arbitrary units, and high 
expressing cells had a total fluorescence greater than 70 units. 

Toxicity assays 

Stock solutions of 50 mM juglone (Calbiochem) and 20 mM 
thimerosal (Enzo) were freshly dissolved in DMSO or water, 



respectively, prior to addition to molten NGM. Sodium arsenite 
(Ricca) was maintained in aqueous solution at 0.5% w/v and 
stored at room temperature. Plates were freshly made approxi- 
mately 24 hours before use and seeded with concentrated OP50 
the night before being used. 

To assess longevity, age matched young adult animals were 
transferred to 100 mm NGM plates containing either drug or 
control and assayed over time. Animals which escaped the plates 
were excluded from the analysis. At least four replicates of n = 40 
animals per genotype per stress were tested; final samples sizes 
reported in Table S3. Animals were stored at 20° except during 
time points. Animals were scored as dead if they did not respond to 
repeated light prodding. Percentages alive for each genotype were 
determined by averaging the fraction alive per plate at each time 
point and plotting graphically. 

For fluorescence toxicity studies using the nlg-1 reporters, L4 
stage animals were exposed to NGM plates supplemented with 
drug or control overnight for 14 hours; animals were allowed 2— 
4 hours recovery time before imaging. Concentrations were 
chosen that did not result in animal death after 14 hours. For 
juglone imaging, control plates were supplemented with an equal 
volume of DMSO, as the juglone was diluted to 50 mM in DMSO 
prior to addition to the test plates. 

Statistical analysis 

A Student's t test was used to determine significance when 
comparing fluorescence of nlg-1 reporters in different conditions, 
unless otherwise specified. Log rank tests with a Bonferroni 
correction were calculated by JMP Pro version 10.0 software and 
were used to determine significant differences between genotypes 
for toxicity studies; differences between genotypes are reported in 
Table S3. 

RNA sequencing 

Total RNA was isolated from approximately 10,000 mixed 
stage animals for wild type, wdr-23(tml817) mutants and wdr- 
23;mh225 using Stat60 (Tel-test B, Texas). Transcriptome 
libraries were prepared using TruSeq RNA sample preparation 
kit (Illumina) according to manufacturer's TruSeq protocol as 
previously described [20]. Libraries were amplified by PCR and 
quality and quantity of libraries were evaluated on BioAnalyzer 
2100 (Agilent). Sequencing was performed on HiSeq2000 
(Illumina). Sequencing reads were aligned to the C. elegansgenome 
(release WS210) using TopHat [68]. Gene models were down- 
loaded from ENSEMBL and quantified using Cufflinks. Differen- 
tially expressed genes at false discovery rate (FDR) of 0.05 were 
identified using the Cuffdiff module of the Cufflinks package. 

Electrophoretic mobility shift assays 

Full-length SKN-1 a cDNA was cloned into a pBS backbone 
driven by the T7 promoter and expressed using TnT T7 Quick 
Coupled Transcription/Translation System for DNA (Promega). 
EMSA was performed using the LightShift Chemiluminescence 
EMSA Kit (Pierce) according to manufacturer's protocols. 
Complementary 5' biotinylated oligonucleotides containing the 
pnlg-1 SKN-1 binding site were self-annealed and incubated with 
1 ul of SKN-1 lysate, 1 ng Poly (dldC) and 5 mM MgCl 2 for 
20 minutes at room temperature. Samples were separated on a 
5% native polyacrylamide gel and blotted on Biodyne B nylon 
membranes (Pierce). 

Pnlg-1 probes: 

oTS300: 5' biotin-gttgccccccaaatgATGACATTacctcttttcctccc 

3' 
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oTS310: 5' biotin-gggaggaaaagaggtAATGTCATcatttgglggg- 
caac 3' 

Supporting Information 

Figure SI nlg-1 promoter alignment for Caenorhabditis species. 
***Indicates transcriptional start aligned for all four species. 
Consensus sequences (WWTDTCAT) were identified by RSAT 
and indicated for C. ekgans by bold, underlined sequences. 
Conserved SKN-1 consensus sequence boxed. 
(TIF) 

Figure S2 Thimerosal increases intestinal SKN-1 ::GFP. Intesti- 
nal fluorescence of SKN-lb/c::GFP in control animals (top) 
compared to animals treated with 35 uM thimerosal for 14 hours. 
Animals were imaged at L3/L4 stage to limit background 
fluorescence. 
(TIF) 

Figure S3 Quantification of Pnlg-lr.gfp (vjls47) in the ventral cord 
of animals in different genetic backgrounds. Samples sizes shown. 
*/><0.05, **/><0.01. 
(TIF) 

Figure S4 Survival in response to sodium arsenite. Survival 
curves of indicated strains on 2.5 mM sodium arsenite. Error bars 
represent ±sem. 
(TIF) 

Figure S5 Survival in response to thimerosal. Survival curves of 
indicated strains on 50 uM thimerosal. Error bars represent 
±sem. 
(TIF) 

Figure S6 Survival of skn-l(gf) on 200 u,M juglone. Survival 
curves of indicated strains on 200 uM juglone. Error bars 
represent ±sem. 
(TIF) 

Figure S7 Survival of Pnlg-lr.gfp on 200 u,M juglone. Survival 
curves of indicated strains on 200 uM juglone. Error bars 
represent ±sem. 
(TIF) 

Table SI RNA sequencing of wdr-23 mutants and WDR- 
23 neuronal rescue. RNAseq list summarizing the 2,285 genes 

References 

1. Finkcl T, Holbrook NJ (2000) Oxidants, oxidative stress and the biology of 
ageing. Nature 408: 239-247. 

2. Halliwell B (2006) Oxidative stress and neurodegeneration: where are we now? 
J Ncurochcm 97: 1634-1658. 

3. Abou-Slciman PM, Muqit MM, Wood NW (2006) Expanding insights of 
mitochondrial dysfunction in Parkinson's disease. Nat Rev Neurosci 7: 207-219. 

4. Zana M, Janka Z, Kalman J (2007) Oxidative stress: a bridge between Down's 
syndrome and Alzheimer's disease. Neurobiol Aging 28: 648—676. 

5. Itoh K, Chiba T, Takahashi S, Ishii T, Igarashi K, et al. (1997) An Nrf2/small 
Maf hetcrodimer mediates the induction of phase II detoxifying enzyme genes 
through antioxidant response elements. Biochem Biophys Res Gommun 236: 
313-322. 

6. Kobayashi M, Yamamoto M (2006) Nrf2-Keapl regulation of cellular defense 
mechanisms against eleetrophiles and reactive oxygen species. Adv Enzyme 
Regul 46: 113-140. 

7. Sykiotis GP, Bohmann D (2010) Stress-activated cap'n'collar transcription 
factors in aging and human disease. Sei Signal 3: re3. 

8. Lee JM, LiJ, Johnson DA, Stein TD, Kraft AD, et al. (2005) Nrf2, a multi-organ 
protector? FASEBJ 19: 1061-1066. 

9. Itoh K, Wakabayashi N, Katoh Y, Ishii T, Igarashi K, ct al. (1999) Keapl 
represses nuclear activation of antioxidant responsive elements by Nrf2 through 
binding to the amino-terminal Neh2 domain. Genes Dev 13: 76-86. 

10. Kraft AD, Johnson DA, Johnson JA (2004) Nuclear factor E2-rclatcd factor 2- 
dependent antioxidant response element activation by tert-butylhydroquinonc 
and sulforaphanc occurring preferentially in astrocytes conditions neurons 
against oxidative insult. J Neurosci 24: 1101-1112. 



up-regulated in wdr-23 mutants. Neuron rescue denotes wdr-23 
mutants expressing an integrated wdr-23a driven by the snb-1 
promoter (vjh225). Genes highlighted in green are statically 
different between wdr-23 mutants and WDR-23 rescue. Expression 
was identified for genes which are significandy different between 
wdr-23 and WDR-23 rescue using WormMart. WWTDTCAT 
denotes SKN-1 binding sites identified within 1000 bp of gene 
start. Genes located proximally to ChlP-seq sites are identified by 
stage [34]; sites were restricted to 2 kb upstream or 500 bp 
downstream of start sites. Genes identified by microarrays are 
identified in Ref column as follows: skn-1 RNAi 1 , exposure to 
arsenite 2 , exposure to <-butyl hydrogen peroxide 3 [42] or 
hyperbaric oxygen 4 [49]. 

(xisx) 

Table S2 Punctal analysis of NLG-1 -GFP fluorescence. Raw 
values ± sem reported for peak, cord, width, and interpunctal 
interval (IPI). 
(PDF) 

Table S3 Statistical analysis of toxicity assays. Samples sizes, % 2 
values for log rank tests, and significant differences between 
samples reported. Statistically significant differences were calcu- 
lated using the log rank test with a Bonferroni correction to 
establish the thresholds of significance. 
(PDF) 

Acknowledgments 

We thank J. Chan and H. Wang for advice and critical review of the 
manuscript, and I. Antoshechkin, in association with the Millard and 
Muriel Jacobs Genetics and Genomics Laboratory at California Institute of 
Technology, for RNAseq analysis and help identifying overlapping ChlP- 
seq peaks. We thank J. Shaw for her technical assistance, and Q. Ch'ng 
and J. Dittman for fluorescence analysis software. We thank S. Curran for 
skn-1 (gf) alleles, L. Dreier for the mito-GFP strain, and A. Van der Bliek for 
the INVOM-RFP construct. 

Author Contributions 

Conceived and designed the experiments: TAS DS. Performed the 
experiments: TAS OE DS. Analyzed the data: TAS DS. Contributed 
reagents/materials/analysis tools: TAS OE JAK DS. Wrote the paper: 
TAS DS. 



11. Lee JM, Shih AY, Murphy TH, Johnson JA (2003) NF-E2-related factor-2 
mediates neuroprotection against mitochondrial complex I inhibitors and 
increased concentrations of intracellular calcium in primary cortical neurons. 
J Biol Chem 278: 37948-37956. 

12. LccJM, Calkins MJ, Chan K, Kan YW, Johnson JA (2003) Identification of the 
NF-E2-related factor-2-dependcnt genes conferring protection against oxidative 
stress in primary cortical astrocytes using oligonucleotide microarray analysis. 
J Biol Chem 278: 12029-12038. 

13. An JH, Blackwcll TK (2003) SKN-1 links C. elcgans mcscndodcrmal 
specification to a conserved oxidative stress response. Genes Dev 17: 1882-1893. 

14. Bishop NA, Guarente L (2007) Two neurons mediate diet-restriction-induced 
longevity in C. clegans. Nature 447: 545-549. 

15. Vanduyn N, Settivari R, Wong G, Nass R (2010) SKN-1/Nrf2 inhibits 
dopamine neuron degeneration in a Caenorhabditis clegans model of 
methylmercury toxicity. Toxicol Sci 118: 613—624. 

16. Settivari R, Vanduyn N, LcvoraJ, Nass R (2013) The Nrf2/ SKN-1 -dependent 
glutathione S-transferase pi homologue GST-1 inhibits dopamine neuron 
degeneration in a Caenorhabditis elcgans model of manganism. Ncurotoxicol- 
ogy 38C: 51-60. 

17. VanDuyn N, Settivari R, LcVora J, Zhou S, Unrine J, ct al. (2013) The metal 
transporter SMF'-S/DMT-l mediates aluminum-induced dopamine neuron 
degeneration. J Neurochem 124: 147-157. 

18. Benedetto A, Au C, Avila DS, Milatovic D, Aschner M (2010) Extracellular 
dopamine potentiates mn-induccd oxidative stress, lifespan reduction, and 
dopaminergic neurodegeneration in a BLI-3-depcndent manner in Caenorhab- 
ditis elegans. PLoS Genet 6: el 00 1084. 



PLOS Genetics | www.plosgenetics.org 



13 



January 2014 | Volume 10 | Issue 1 | e1004100 



Stress Regulation of Neuroligin 



19. Choc KP, Przybysz AJ, Strange K (2009) The WD40 repeat protein WDR-23 
functions with the CUL4/DDB1 ubiquitin ligase to regulate nuclear abundance 
and activity of SKN-1 in Caenorhabditis elegans. Mol Cell Biol 29: 2704-2715. 

20. Staab TA, Griffen TC, Corcoran C, Evgrafov O, Knowlcs JA, et al. (2013) The 
Conserved SKN-1 /Nrl2 Stress Response Pathway Regulates Synaptic Function 
in Caenorhabditis elegans. PLoS Genet 9: el 003354. 

21. Hascgawa K, Miwa J (2010) Genetic and cellular characterization of 
Caenorhabditis elegans mutants abnormal in the regulation of" many phase II 
enzymes. PloS one 5: c 1 1 1 94. 

22. wormbasc.org Wormbasc web site. Release \vs236 ed. 

23. Blumenthal T (2012) Trans-splicing and opcrons in C. elegans. WormBook: 1- 
11. 

24. Melcndcz A, Talloczy Z, Seaman M, Eskclinen EL, Hall DH, et al. (2003) 
Autophagy genes are essential for dauer development and life-span extension in 
C. elegans. Science 301: 1387-1391. 

25. Takacs-Vellai K, Vellai T, Puoti A, Passannante M, Wicky C, et al. (2005) 
Inactivation of the autophagy gene bec-1 triggers apoptotic cell death in C. 
elegans. Curr Biol 15: 1513—1517. 

26. Hall DH, Hcdgecock EM (1991) Kinesin-related gene une-104 is required for 
axonal transport of synaptic vesicles in C. elegans. Cell 65: 837—847. 

27. Labrousse AM, Zappaterra MD, Rube DA, van der Bliek AM (1999) C. elegans 
dynamin-related protein DRP-1 controls severing of the mitochondrial outer 
membrane. Mol Cell 4: 815-826. 

28. Dennis G, Jr., Sherman BT, Hosack DA, Yang J, Gao W, et al. (2003) DAVID: 
Database for Annotation, Visualization, and Integrated Discovery. Genome Biol 
4: P3. 

29. Iwasaki K, Toyonaga R (2000) The Rab3 GDP/GTP exchange factor homolog 
AEX-3 has a dual function in synaptic transmission. EMBO J 19: 4806-4816. 

30. Hunter JW, Mullen GP, McManus JR, HeatherlyJM, Duke A, ct al. (2010) 
Neuroligin-dcficicnt mutants of C. elegans have sensory processing deficits and 
are hypersensitive to oxidative stress and mercury toxicity. Dis Model Mech 3: 
366-376. 

3 1 . Kohn RE, Duerr JS, McManus JR, Duke A, Rakow TL, et al. (2000) Expression 
of multiple UNC-13 proteins in the Caenorhabditis elegans nervous system. Mol 
Biol Cell 11: 3441-3452. 

32. Charlie NK, Thomurc AM, Schadc MA, Miller KG (2006) The Dunce cAMP 
phosphodiesterase PDE-4 negatively regulates G alpha(s)-depcndent and G 
alpha(s) -independent cAMP pools in the Caenorhabditis elegans synaptic 
signaling network. Genetics 173: 111-130. 

33. Thomas-Chollicr M, Sand O, Turatsinze JV, Janky R, Defrance M, et al. (2008) 
RSAT: regulatory sequence analysis tools. Nucleic Acids Res 36: Wl 19-127. 

34. Niu W, Lu ZJ, Zhong M, Sarov M, Murray JI, et al. (201 1) Diverse transcription 
factor binding features revealed by genome-wide ChlP-seq in C. elegans. 
Genome Res 21: 245-254. 

35. Hu Z, Horn S, Kudze T, Tong XJ, Choi S, et al. (2012) Neurexin and neuroligin 
mediate retrograde synaptic inhibition in C. elegans. Science 337: 980-984. 

36. Calahorro F, Alcjandre E, Ruiz-Rubio M (2009) Osmotic avoidance in 
Caenorhabditis elegans: synaptic function of two genes, orthologues of human 
NRXN1 and NLGN1, as candidates for autism. J Vis Exp (34): 1616. 

37. Feinberg EH, Vanhoven MK, Bendesky A, Wang G, Fetter RD, et al. (2008) 
GFP Reconstitution Across Synaptic Partners (GRASP) defines cell contacts and 
synapses in living nervous systems. Neuron 57: 353-363. 

38. Pack J, Lo JY, Narasimhan SD, Nguyen TN, Glover-Cutter K, et al. (2012) 
Mitochondrial SKN-1 /Nrf Mediates a Conserved Starvation Response. Cell 
Mctab 16: 526-537. 

39. Blackwell TK, Bowerman B, Priess JR, Weintraub H (1994) Formation of a 
monomeric DNA binding domain by Skn-1 bZIP and homeodomain elements. 
Science 266: 621-628. 

40. Sieburth D, Madison JM, Kaplan JM (2007) PKC-1 regulates secretion of 
neuropeptides. Nat Neurosci 10: 49-57. 

41. Dittman JS, Kaplan JM (2006) Factors regulating the abundance and 
localization of synaptobrevin in the plasma membrane. Proc Natl Acad 
Sci USA 103: 11399-11404. 

42. Oliveira RP, Porter Abate J, Dilks K, LandisJ, Ashraf J, ct al. (2009) Condition- 
adapted stress and longevity gene regulation by Caenorhabditis elegans SKN-1/ 
Nrf. Aging Cell 8: 524-541. 

43. Przybysz AJ, Choc KP, Roberts LJ, Strange K (2009) Increased age reduces 
DAF-16 and SKN-1 signaling and the hormetic response of Caenorhabditis 
elegans to the xenobiotie juglone. Mech Ageing Dev 130: 357-369. 



44. Tullet JM, Hertweck M, An JH, Baker J, Hwang JY, et al. (2008) Direct 
inhibition of the longevity-promoting factor SKN- 1 by insulin-like signaling in 
C. elegans. Cell 132: 1025-1038. 

45. Inoue H, Hisamoto N, An JH, Oliveira RP, Nishida E, et al. (2005) The C. 
elegans p38 MAPK pathway regulates nuclear localization of the transcription 
factor SKN-1 in oxidative stress response. Genes Dev 19: 2278-2283. 

46. Zubovych IO, Straud S, Roth MG (2010) Mitochondrial dysfunction confers 
resistance to multiple drugs in Caenorhabditis elegans. Mol Biol Cell 21: 956- 
969. 

47. Narendra D, Tanaka A, Suen DF, Youle RJ (2008) Parkin is recruited selectively 
to impaired mitochondria and promotes their autophagy. J Cell Biol 183: 795— 
803. 

48. Youle RJ, van der Bliek AM (2012) Mitochondrial fission, fusion, and stress. 
Science 337: 1062-1065. 

49. Park SK, Tedcsco PM, Johnson TE (2009) Oxidative stress and longevity in 
Caenorhabditis elegans as mediated by SKN-1. Aging cell 8: 258—269. 

50. Luo L, O'Lcary DD (2005) Axon retraction and degeneration in development 
and disease. Annual review of neuroscienee 28: 127—156. 

51. Feng Z, Li L, Ng PY, Porter AG (2002) Neuronal differentiation and protection 
from nitric oxide-induced apoptosis require c-Jun-dcpendcnt expression of 
NCAM140. Mol Cell Biol 22: 5357-5366. 

52. Klementicv B, Novikova T, Korshunova I, Bcrezin V, Bock E (2008) The 
NCAM -derived P2 peptide facilitates recovery of cognitive and motor function 
and ameliorates neuropathology following traumatic brain injury. EurJ Neurosci 
27: 2885-2896. 

53. Scapagnini G, Vasto S, Abraham NG, Caruso C, Zella D, et al. (2011) 
Modulation of Nrf2/ARE pathway by food polyphenols: a nutritional 
neuroprotective strategy for cognitive and neurodegenerative disorders. 
Molecular neurobiology 44: 192-201. 

54. Chakraborty S, Aschner M (2012) Altered manganese homeostasis: implications 
for BLI-3-dependent dopaminergic ncurodegeneration and SKN-1 protection in 
C. elegans. Journal of trace elements in medicine and biology : organ of the 
Society for Minerals and Trace Elements 26: 183—187. 

55. Joshi G, Johnson JA (2012) The Nrf2-ARE pathway: a valuable therapeutic 
target for the treatment of neurodegenerative diseases. Recent patents on CNS 
drug discovery 7: 218-229. 

56. Bertholet AM, Millet AM, Guillermin O, Daloyau M, Davezac N, ct al. (2013) 
OPA1 loss of function affects in vitro neuronal maturation. Brain : a journal of 
neurology' 136: 1518-1533. 

57. DurieuxJ, Wolff S, Dillin A (2011) The cell-non-autonomous nature of electron 
transport chain-mediated longevity. Cell 144: 79—91. 

58. Sudhof TC (2008) Ncuroligins and neurcxins link synaptic function to cognitive 
disease. Nature 455: 903-911. 

59. Bottos A, Rissonc A, Bussolino F, Arese M (2011) Neurcxins and ncuroligins: 
synapses look out of the nervous system. Cell Mol Life Sci 68: 2655-2666. 

60. Zhang C, MilunskyJM, Newton S, Ko J, Zhao G, et al. (2009) A neuroligin-4 
missense mutation associated with autism impairs neuroligin-4 folding and 
endoplasmic reticulum export. J Neurosci 29: 10843-10854. 

61. Comoletti D, De Jaco A, Jennings LL, Flynn RE, Gaietta G, et al. (2004) The 
Arg45 1 Cys-neuroligin-3 mutation associated with autism reveals a defect in 
protein processing. J Neurosci 24: 4889—4893. 

62. Chih B, Engclman H, Schciffclc P (2005) Control of excitatory and inhibitory 
synapse formation by ncuroligins. Science 307: 1324—1328. 

63. Ghezzo A, Visconti P, Abruzzo PM, Bolotta A, Ferrcri C, et al. (2013) Oxidative 
Stress and Erythrocyte Membrane Alterations in Children with Autism: 
Correlation with Clinical Features. PLoS One 8: e66418. 

64. Gorrindo P, Lane CJ, Lee EB, McLaughlin B, Levitt P (2013) Enrichment of 
elevated plasma F2t-isoprostanc levels in individuals with autism who are 
stratified by presence of gastrointestinal dysfunction. PLoS One 8: c68444. 

65. Rossignol DA, Fryc RE (2012) Mitochondrial dysfunction in autism spectrum 
disorders: a systematic review and meta-analysis. Mol Psychiatry 17: 290-314. 

66. Mello CC, KramerJM, Stinchcomb D, Ambros V (1991) Efficient gene transfer 
in C. elegans: cxtrachromosomal maintenance and integration of transforming 
sequences. EMBO J 10: 3959-3970. 

67. Ch'ng Q> Sieburth D, Kaplan JM (2008) Profiling synaptic proteins identifies 
regulators of insulin secretion and lifespan. PLoS Genet 4: el000283. 

68. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, et al. (2010) 
Transcript assembly and quantification by RNA-Seq reveals unannotated 
transcripts and isoform switching during cell differentiation. Nat Biotechnol 28: 
511-515. 



PLOS Genetics | www.plosgenetics.org 



14 



January 2014 | Volume 10 | Issue 1 | e1004100 



